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Abstract 



Earlier works in engineering, partly experimental, partly computational have revealed that asymptotically, 
when the excitation is a white noise, plastic deformation and total deformation for an elasto-perfectly-plastic 
oscillator have a variance which increases linearly with time with the same coefficient. In this work, we prove 
this result and we characterize the corresponding drift coefficient. Our study relies on a stochastic variational 
inequality governing the evolution between the velocity of the oscillator and the non-linear restoring force. 
We then define long cycles behavior of the Markov process solution of the stochastic variational inequality 
which is the key concept to obtain the result. An important question in engineering is to compute this 
coefficient. Also, we provide numerical simulations which show successful agreement with our theoretical 
prediction and previous empirical studies made by engineers. 



Resume 

Le comportement de la deformation plastique pour un oscillateur elastique-parfaitement- 
plastique excite par un bruit blanc. Des resultats experimentaux en sciences de l'ingenieur ont montre 
que, pour un oscillateur elasto-plastique-parfait excite par un bruit blanc, la deformation plastique et la de- 
formation totale ont une variance, qui asymptotiquement, croit lineairement avec le temps avec le meme co- 
efficient. Dans ce travail, nous prouvons ce resultat et nous caracterisons le coefficient de derive. Notre etude 
repose sur une inequation variationnelle stochastique gouvernant revolution entre la vitesse de l'oscillateur 
et la force de rappel non-lineaire. Nous definissons alors le comportement en cycles longs du processus de 
Markov solution de Pinequation variationnelle stochastique qui est le concept cle pour obtenir le resultat. 
Une question importante en sciences de l'ingenieur est de calculer ce coefficient. Les resultats numeriques 
confirment avec succes notre prediction theorique et les etudes empiriques faites par les ingenieurs. 



Version francaise abregee 

Dans cet article, nous etudions la variance de l'oscillateur elastique-parfaitement-plastiquc (EPP) excite 
par un bruit blanc. La dynamique de l'oscillateur s'exprime a l'aide d'une equation a memoire. (voir 
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(H}-©). A.Bensoussan and J.Turi ont montre que la relation entre la vitesse et la composante elastique 
satisfait une inequation variationnelle stochastique (voir (| 1 SVZ[) '). Dans ce cadre, nous introduisons les cycles 
long independants (definis plus loin) et nous justifions qu'ils permettent de caracteriser la variance de la 
deformation totale et de la deformation plastique (voir Q). 



1. Introduction 

In civil engineering, an elasto-perfectly-plastic (EPP) oscillator with one single degree of freedom is 
employed to estimate prediction of failure of mechanical structures under random vibrations. This elasto- 
plastic oscillator consists in a one dimensional model simple and representative of the elasto-plastic behavior 
of a class of structure dominated by their first mode of vibration. The main difficulty to study these systems 
comes from a frequent occurrence of plastic phases on small intervals of time. A plastic deformation is 
produced when the stress of the structure crosses over an elastic limit. The dynamics of the EPP-oscillator 
has memory, so it has been formulated in the engineering literature as a process with hysteresis x(t), which 
stands for the displacement of the oscillator. We study the problem 

x + c ± + Ft = w (1) 

with initial conditions of displacement and velocity x(0) — 0, i(0) = 0. Here Co > is the viscous damping 
coefficient, k > the stiffness, w is a Wiener process; F t := F(x(s),0 < s < t) is a non-linear functional 
which depends on the entire trajectory {x(s), < s < t} up to time t. The plastic deformation denoted by 
A(t) at time t can be recovered from the pair (x(t), F t ) by the following relationship: 

( kY if x(t)=Y + A(t), 

F t = \ k(x(t)-A(t)) if x(t)e]-Y + A(t),Y + A(t)[, (2) 
[ -kY if x(t) = -Y + A(t). 

where A(t) = y(s)l{|F t |=fcy}ds and Y is an elasto-plastic bound. Karnopp & Scharton [j| proposed a 
separation between elastic states and plastic states. They introduced a fictitious variable z{t) := x(t) — A(t) 
and noticed the simple fact that between two plastic phases z(t) behaves like a linear oscillator. In addition, 
other works made by engineers [y], partly experimental, partly computational, have revealed that total 
deformation has a variance which increases linearly with time: 

lim^^W. (3) 

t— s-oo t 



where a 2 is a positive real number. 



Recently, the right mathematical framework of stochastic variational inequalities (SVI) modeling an 
elasto-plastic oscillator with noise (presented below) has been introduced by the first author with J. Turi in 
[3]. The inequality governs the relationship between the velocity y(t) and the elastic component z{i): 

dy(t) = -(c y(t) + kz{t))dt + dw(t), (dz(t)-y(t)dt)(<t>-z(t))>0, V|^| < Y, \z(t)\ < Y. {SVI) 

The plastic deformation A(t) can be recovered by J* y(s)l^ z ^ s ^ = Y}ds. 

Throughout the paper, the objective is to prove ([3]) and to provide a exact characterization of §3§ based 
on (|tSVT[) . Now, we introduce long cycles in view to the fact that we are interested in identifying independent 
sequences in the trajectory. 
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1.1. Long cycles 

Denote tq :— mi{t > 0, y(t) = and \z(t)\ = Y} and s := sign(z(ro)) which labels the first 
boundary hit by the process (y(t), z(t)). Define 6q := inf{< > To, y(t) = and z(t) = —sY}. In a 
recurrent manner, knowing 9 n , we can define for n > 

r n+1 := M{t > 9 n , y(t) = Q and z(t) = sY}, 
n+ i := inf{t > r n+ i, y(t) = and z(t) = -sY}. 

Now, according to the previous setting we can define the n-th long cycle (resp. first part of the cycle, 
second part of the cycle) as the piece of trajectory delimited by the interval [r n , T n +i), (resp. [r n , and 
[0 ra +i, t„ + i)). Indeed, at the instants {T ni n > 1}, the process z(t)) is in the same state at the instant 

tq. In addition, there are two types of cycles depending on the sign of s. The set of stopping times {r„, n > 0} 
represents the times of occurence of long cycles. Let us remark that the plastic deformation and the total 
deformation are the same on a long cycle since J^ 1 y(t)dt = J^ 1 y(t)l{\ z ( s )\ = Y}dt + J^ 1 ytyl^g^^yydt and 
that J^ 1 y(t)l{| 2 ( s )|<y}dt = z(t\) — z(tq) = 0. As main result, we have obtained the following theorem. 

Theorem 1.1 (Characterization of the variance related to the plastic/total deformation). In the 

previously defined context, we have shown 

*nm = e(qm*) 2 

f^oo t E(n - t ) 

Our proofs are based on solving nonlocal partial differential equations related to long cycles. Simpler formula 
will be given below at equation (|13D . 

2. The issue of long cycles and plastic deformations 

Let us introduce notations. 

Notation 1. D := R x (-Y,+Y), D+ := (0, oo) x {Y}, D~ := (-oo,0) x {-Y} ; and the differential 
operators A( := -±( yy + (c Q y + kz)( y -y( z , B + ( := -\Cyy + (coy + kY)( y , B-( := -\C, yy + {c y-kY)C, y . 
where Q is a regular function on D. 

Let / be a bounded measurable function, we want to solve 

Av+=f(y,z) in D, B+v+ = f(y,Y) in D+, B_v+ = f(y,-Y) in D~~ (P„+) 
with the nonlocal boundary conditions v + (y, Y) continuous and u + ((Y, — Y) = 0, and 

Av-=f(y,z) in D, B+v~ = f(y,Y) in D+, B-v~ = f(y,-Y) in D~ (P v - ) 

with the nonlocal boundary conditions v~(0 + , Y) = and v~(y, —Y) continuous. The functions v + (y, z\ f) 
and v~(y,z;f) are called long cycles. In addition, we define 7r + (y, z) and ir~(y,z) by 

Att + =0 in D, n + (y,Y) = l in D + , vr + (y, -Y) = in IT (P^+) 

and 

Air' =0 in D, 7r _ (y,Y)=0 in D + , n~(y,-Y) = l in £T (P,-) 

Note that ir + (y, z) + ir~ (y, z) — 1, so the existence and uniqueness of a bounded solution \P^+ 1 and \P^- 1 
is clear. 

Proposition 2.1. IVe /iai/e the properties 7r~(0 + , — Y) > and 7r + ((Y, Y) > 0. 
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Proof. We check only the first property. We consider the elastic process (y yz (t), z yz (t)): 

z v Jt) = e~ ^~{z cos (ut) + — (y + —z) sin (uit)} + — f e~^ t_s ' sin (u(t — s))dw(s), 

oj 2 uj J 

rt 

Vyz(t) = —-?rZyz(t) + e 2 {— wzsin (ujt) + (y + —z) cos (u)t)} + e~~ 2 ( t-a ) cos (u(t — s))dw(s) 
2 2 7n 



where u> := c ° (we assume 4fc > Cq). Note that the condition 4fc > Cq is needed so that (y yz (t), z yz (t)) 
have real valued solutions. Set r yz := inf{£ > 0, > Y} then we have the probabilistic interpretation 

TT+(y,z) = P(z yz (T yz ) = Y), 7T~(y,z) = f(z yz (T yz ) = -Y). We can state 

7r _ (y, z) —> 1 as i/^-oo, z £ [—Y,Y]. (5) 

Indeed, Vi with < t < — we have z yz (t) — > — oo, as y — > — oo a.s. Therefore Vt with < i < — , 
a.s. z yz (t) < —Y for y sufficiently large. Hence, a.s. r yz < t for y sufficiently large. Therefore a.s. 
lim sup J) _ > _ 00 r yz < t. Since t is arbitrary, necessarily a.s. r yz — > 0, as y — > —oo which implies ((S]). Moreover 
the function 7r + (y, Y) cannot have a minimum or a maximum at any finite y < 0. It is then monotone and 
since ir~(— oo,Y) = 1, it is monotone decreasing. It follows that ir~(0~,Y) < 1. It cannot be 0. Otherwise, 
7r~(y, Y) is continuous at y = 0, and (0, Y) is a point of minimum of 7r~ (y, z). Since for y < 0, 7T~ (y, Y) < 
from the equation of 7r _ we get limsup^Q Tr~ y (y, Y) < which is not possible since (0, Y) is a minimum. 

We next define r)(y,z) by 
1 
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i] yy + (c y + kz)ri y - yi] z = f(y,z) in D, ??(y,Y)=0 in D + , r?(y, — Y") = in D (P, 



with the local boundary conditions r)(0 + ,Y) — and r;(0 , — Y) = 0. For / bounded (P n ) has a unique 
bounded solution. Then define (p+(y, /) by solving 

-~ip+, y y + (c y + kY)ip + , y =f(y,Y), y > 0, p+(0;/) = 0. (6) 
We check easily the formula 

poo r£+y 

V+(y;f) = 2 d£exp(-(co£ 2 + 2£YO) / f(C, Y)exp(-2c £(C - 0)dC, V > 0. (7) 
Jo J£ 

We consider ip+iy, z; /) defined by 

-^ip +t y y + (c y + kz)ip +t y-yip +tZ = in D, 

^+(y,Y) = v + {y:f) in D+, (P^) 
1>+(V,-Y) = in D~ 



and similarly ip—(y; /), ip—(y, z\ f) defined by 
1 

~ 2 



^>-, M + (c y-»>_ i3 , = /(y,-Y), y<0, y>_(0;/) = 0. (8) 



which leads to 



/•oo /•- s 

<p-(y;f) = 2 d£cxp(-( Co £ 2 -2fcYO) / /(C; -V) exp(-2c £(C - OMC, V < 0. 
Jo Jy-J 

and 

+ (coy + kz)ip-, y - yip-, z = in D, 



V>-(y,Y) - in D+, (P^J 

^-(y,-^) = <p-(v''f) in D ~ 



We can state 



Proposition 2.2. The solution of (K+J is given by 



v + (y,z;f) =r)(y,z;f) + ip + (y,z;f) + tp-(y,z;f) (9) 
, f?(Q-,y; /) + ^+(o-,y ; /) + il>-(0-,Y; f) 



7T-(0-,F) 



-7r + (y,z). 



and the solution of jP v - D is given by 



v (y,z;f)=r)(y,z;f) + ip + (y,z;f) + ip-(y,z;f) (10) 
77(0+, -Y; /) + V+(0+, — Y; /) + V-(0+, -Y; /) 



7T+(0-,Y) 

Proof. Direct checking. 



-7T (y,s). 



3. Complete cycle 

First, let us check that E[x(t)] = and then a 2 (x(t)) = K[x 2 (t)}. Indeed, by symmetry of the inequality 
(16*yZ[) and by the choice of initial conditions y(0) = 0, z(0) = 0, the processes (y(t),z(t)) and —(y(t),z(t)) 

have same law. Then, E[a;(t)] = E f Q y(s)ds = 0. In addition, (y(ri), z(t\)) is equal to (0, — Y) or 
(0, Y) with probability 1/2 for both, therefore with no loss of generality we can suppose that (y(0),z(0)) = 
(0, -Y) or (0,Y) with probability 1/2 for both. Let us treat the case y(0) = 0,2(0) = Y. So, t = 
and 9i = inf{/j > 0, z(t) = -Y,y(t) = 0}. We can assert that E6>i = w+(0,Y;l) hence 6 1 < oo 
a.s.. Next we define t x = inf{t > 9 X , z(t) = Y,y(t) = 0} then Er x = u+(0,Y;l) + tr (0, - Y; 1) = 
2u + (0,Y;l). At time T\ the state of the system is again (0, Y). So the sequence {r n ,n > 0} is such 
that r n < t„ + i and in the interval (r„,r„ + i) we have a cycle identical to (0,ri). Consider the variable 
J T1 y{t)dt. We have E /J" 1 y(t)dt = E tf 1 y(t)dt + E j/(i)tft = w+ (0, Y; y) + tr (0, - Y; y). However if / 
is antisymmetric f(—y,—z) = —f(y,z) we have w + (0,Y; /) = — v~(0, — Y; /), therefore E JJ" 1 y(t)dt = 0. 

Hence, E (/ T " y(i)dt) 2 = E (^d J^ +1 l/(«)di) * = nE (J^ y(i)dt) 2 . Then 



E(J^y(t)dt) 2 _E(£y(t)dtY 



Et„ En 

Let T > and iVr with tjv t < T < tjv t+1 , iVy = if Ti > T, ro = and then calculations lead to 
the following 

E(j Q T ^y(t)dt) 2 _E(j^y(t)dt) 2 



Next, we can justify that 




Etjv t+ i Eti 



1 



1~N T + 1 



y(t)dtj = T lm-E^ y(t)dtj (11) 
and that we have a lower bound and an upper bound for pip , that is 

'"y' < ±E ( 2 < f i + ^ E(J °' i "" d ' )2 

Therefore 



1 / f T ,s V E(f n Tl y(t)dt) 2 , s 

T lim -E / y(t)dt\ = UoVU ) , (12) 



Moreover, we can simplify (|12[) . indeed we have 



EU y(t)dtj =E 
and therefore, we can justify that 



y(t)dt\ + (J^ y(t)dtj +2J^ y(t)dt I ,/(/),!/ 



p n \ 2 

E( J y(t)dt) =2 

— 0„+ 




y(i)dt -(/ 2/(i)di) 5 
o J Jo 



= 2 




t/(t)dt -(v+(0,Y;y)) 2 



Since Eri = 2w + (0, Y; 1) we obtain 



1 



lim -E| / y(t)dt) = 



E(7 %(i)di) 2 -(v+(0,Y;y)) 2 
v+(0,Y; 1) 



(13) 



4. Numerical evidence in support of our result 

In this section, we provide computational results which confirm our theoretical results. A C code has 
been written to simulate (y(t), z(t)). See Q for the numerical scheme considered to do direct simulation. 
Let T > 0, TV G N and {t n = nSt, < n < N} where St = jj. Then, to compute the left hand side of fl3|, we 
consider MC E N and we generate MC numerical solutions of (|gPT|) {y i (t), < t < T, 1 < i < MC} 



up to the time T . By the law of large numbers, we can approximate ^E ^ y(s)ds^j by X 
TJTc E^? (Eili V l {tn)5t} and also ^E (7 Q T y{s)ds^j by X 2 MC := pjjpE^ (E 
note sx ■— y/Xfrj C — (Xmc) 2 , we also know by the central limit theorem that yE ^ y(s)d 



MC ■ = 

y l {t n )5tj\ Dc- 

2 

* ) e [X M c - 



1.96s 



— - ,Xmc + 1 'f??Z ] with 95% of confidence. Similarly, to compute the right hand side of ((J]), we gen- 
erate MC numerical long cycles. For each trajectory {y l {t),t > 0}, we consider N l c the required num- 



ber of time step to obtain a completed cycle. Denote Smc 

4 



'MC 



and s T :— y/ ^mc 



(tmc) 2 - We also know that ^msl e [ 



1. 963,5 , 1.9633 



r .j 63 I with 95% of confidence. In 
™ c+ -7Wc ™ c ~-7wc 

table [TJ is shown a comparison of the results obtained for Xmc an d where T = 500 , St = 10~ 4 and 
MC = 5000. 



Remark 4.1. From a numerical point of view, the behavior in long cycles is very relevant. Indeed, com- 
puting the left hand side of ^ is much more expensive in time compared to the right hand side (see the 
TMC'-column oftable\7$). 
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Co = 1 , k = 1 


Y 


x MC ,T = mo 




T~MC 


Relative error % 


0.1 


.807±0.03l 


.834±o.o«9 


6 . 61 ±o.u 


3.2 


0.2 


0. 649 ±U.U26 


624 ±0.U4V 


8 . 74 ±U.13 


3.8 


0.3 


.493±o.OiiO 


464 ±U.034 


10.45 ±U - 1B 


5.8 


0.4 


0. 36 i±u.oi4 


0.355 ±UU ^ 


1212 ±u.i8 


1.7 


0.5 


266 ±U.01i 


257±0.019 


13.80 ±0 ^ 


3.3 


0.6 


.195±o.0O8 


O198±() .0i4 


16.15 ±0 - 2ti 


1.5 


0.7 


q 137 ±0.005 


0149 ±o.oii 


lg.84±u.3i 


8 


0.8 


0103 ±U.0U4 


H2 ±0 - 008 


22.80 ±oa9 


8 


0.9 


0. 07 l±o.ooa 


0.086 ±u0,Jti 


26.79 ±U4Y 


15 



Table 1: Monte-Carlo simulations to compare numerical solution of the left and right hand sides of l|4j, T = 500 , St = 10 4 
and MC = 5000. 
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